Optical Polarization of 13 C Nuclei in Diamond through Nitrogen- Vacancy Centers 
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We determine the polarization of the bulk 13 C nuclear spin system in diamond produced by 
interaction with optically oriented nitrogen- vacancy (NV-) defect centers. 13 C nuclei are polarized 
into the higher energy Zeeman state with a bulk-average polarization up to 5.2%, although local 
polarization may be higher. The kinetics of polarization are temperature independent, and occur 
within 5 minutes. Fluctuations in the dipolar field of the NV- center spin bath are identified as 
the mechanism by which nuclear spin transitions are induced near defect centers. Polarization is 
then transported to the bulk material via spin diffusion, which accounts for the observed kinetics 
of polarization. These results indicate control over the nuclear spin bath, a methodology to study 
dynamics of an NV- center ensemble, and application to sensitivity-enhanced NMR. 

PACS numbers: 76.60.-k, 76.30.Mi, 82.56.Na 



In recent years, much attention has been given to the 
negatively charged nitrogen- vacancy (NV-) defect in dia- 
mond. This interest is due to its ground-state triplet, 
long coherence times, large spin polarization via op- 
tical pumping, and optical spin-state readout. Most 
studies have focused on quantum control and the use 
of NV- centers as qubits for quantum informatiorP^, 
while others have focused on applications such as high- 
resolution magnetometry^i. A topic of particular inter- 
est is the interaction of NV- centers with the surrounding 
"spin bath," which contains both nuclear and electronic 
spins!™ 

We have investigated the NV- defect in diamond from 
the perspective of bulk 13 C nuclear polarization in the 
presence of a large, static Zeeman field. Previous work 
has shown microwave-induced dynamic nuclear polariza- 
tion (DNP) of bulk 13 C spins in diamond through param- 
agnetic defect*B"23. These DNP techniques are limited, 
however, by the electron Boltzmann factor at the mag- 
netic field strength and temperature used in the experi- 
ment. Recently, experiments have shown polarization of 
the nitrogen nuclear spin of the NV- center, as well as 
polarization of proximate single 13 C spins by mak ing use 
of the optical polarization of the NV- cente!^^. These 
studies involved polarization and detection of single nu- 
clear spins with optical methods. 

In this Brief Report, we demonstrate an all-optical 
method of polarizing the bulk 13 C (spin^) nuclear spin 
system in diamond at 9.4 Tesla and cryogenic temper- 
atures. We propose a mechanism for the polarization 
which relies upon the fluctuating dipolar field of the NV- 
center spin bath to induce nuclear spin transitions. The 
equilibrium behavior of the polarization is modeled by in- 
voking the spin-temperature concept, where the nuclear 
spin temperature equilibrates with the temperature of 
the dipolar reservoir of the NV- center spin bath. 

The sample used in this study is a 3x3x0. 5mm high 
pressure, high temperature (HPHT) synthetic single- 
crystal diamond purchased from Element-6. After irra- 
diation treatment and annealing, the final concentration 



of NV- centers was 8 ppm (~ 1.4 * 10 18 cm" 3 )Pl 

NMR experiments were carried out in a 9.4 Tesla su- 
perconducting magnet, where spectra were acquired with 
a single-resonance spectrometer at 100.591 MHz. The 
sample was mounted in a homebuilt NMR probe em- 
ploying a flattened split-solenoid coil for optical access. 
A (111) crystal axis was aligned with the magnetic field. 
Temperature control was maintained using an Oxford 
Spectrostat continuous flow helium cryostat with an Ox- 
ford ITC-503S temperature controller. Optical pumping 
was performed using an argon ion laser operating simul- 
taneously at multiple wavelengths ranging from 457. 9nm 
to 514. 5nm. Power output from the laser was 2.5 Watts 
except where laser-power dependence was studied. Opti- 
cal pumping experiments were performed by saturation 
of the 13 C spin transitions via a series of 50 | pulses, each 
followed by a 10 ms dephasing time. After saturation, a 
variable period of optical pumping occurs, then a f pulse 
immediately followed by detection of the free-induction 
decay (FID). Laser irradiation was maintained through- 
out the entire experiment. Non optically-pumped ther- 
mal equilibrium spectra were acquired by cooling the 
sample to the specified lattice temperature, allowing the 
spins to equilibrate for two hours (much longer than the 
nuclear Ti), and applying a f pulse to detect the FID. 
The integrated area of the Fourier transformed FID is 
taken to be proportional to the bulk-average nuclear po- 
larization. Error from random noise was estimated from 
the root-mean-square magnitude and correlation time of 
the noise in the spectrum, and was found to be negligible. 
The error that is apparent in the data with less signal is 
likely due to imperfect baseline and phase corrections. 
Lorentzian curves accurately fit the lineshapes and were 
used to obtain peak shifts and widths. 

In all experiments, optical pumping resulted in 13 C 
NMR signals opposite in sign to thermal equilibrium. 
The magnitude of the polarization increased with lower 
temperature (Fig. [T]) and larger irradiation power (Fig. 

t, with a maximum bulk- average polarization of 5.2%. 
s discussed later, the illuminated volume is less than 
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FIG. 1: (Color online) 13 C polarization kinetics with 2.5W FIG. 2: (Color online) Polarization kinetics as a function of 
irradiation and varied temperature. The theoretical model is laser power at 50K. 
fit to the 5K data. 



the total volume of the sample, so the local polarization 
is significantly greater than 5%. The polarization rate 
appears to be identical in all cases. 

The NMR peak width and frequency shift showed 
a modest temperature dependence, with derived values 
given in Table [TJ It should be noted that significant laser 
heating of the sample may occur, and for temperatures 
20K and lower, the stability of temperature control was 
poor due to the limited cooling capacity of the cryostat. 



Temperature 


Peak Width 


Peak Shift 


5K 


776Hz 


-363Hz 


7.5K 


757Hz 


-315Hz 


10K 


762Hz 


-272Hz 


15K 


731Hz 


-116Hz 


20K 


698Hz 


OHz 



TABLE I: NMR peak width and shift relative to the peak at 
20K, obtained by fitting Lorentzian curves to the data. 

These experimental results give insight into the mech- 
anism of the polarization process. The lack of temper- 
ature dependence of the polarization rate, along with 
the presence of a large Zeeman field and low tempera- 
ture, suggests a process which does not require energy 
exchange with the crystal lattice. We therefore consider 
energy-conserving transitions between 13 C Zeeman and 
NV- center dipolar energy levels. Similar processes have 
been examined as spin-lattice relaxation mechanisms in 
diamo nd an d other materials containing paramagnetic 
centerPESI. 

The total spin Hamiltonian for the system may be writ- 
ten as: 

H = H S +His+ Hi (1) 

Hs refers to terms involving only NV- centers, Hi in- 
volves only nuclear spins, and His represents the in- 
teractions between the two. In this treatment, His is 



a small perturbation that couples the two energy reser- 
voirs independently defined by Us and Hi- We consider 
a simple model that system S is composed of two NV- 
defects. Including pairwise dipolar interactions (between 
defects labeled 1 and 2 ) and retaining only those terms 
that commute with the dominant Zeeman Hamiltonian: 

Hs = hj s (S liZ + S 2iZ ) + T^zf + T 2 {S 2 , Z ? (2) 

+ASi : zS 2 ,z — -^A(Si :+ S 2 - + Si-S 2:+ ) 

where A = j^ 2 7^ryp-[l ~~ 3cos 2 #] quantifies the 
electron-electron dipolar interaction between defects, and 
= 2.88GHz*27rfi.* \ (3 cos 2 74 — l) is the magnitude of 
the spin-spin interaction within a single defeclP^, which 
we have indexed to account for different defect orienta- 
tions. 9 and r specify the orientation and length of the 
vector connecting the two NV- centers and 7 is the an- 
gle of a defect symmetry axis with the external magnetic 
field (0 or 109.5 degrees due to the tetrahedral symmetry 
of the crystal). Diagonalization of Hg gives the energy 
levels for the NV- centers. In order for the perturbation 
His to induce nuclear spin transitions, an accompanying 
transition at the nuclear Zeeman energy must take place 
within the energy levels of Hs- Due to the large Zeeman 
energy (^263 GHz) of NV- centers, transitions between 
NV- center Zeeman levels can never be "on resonance" 
with a nuclear spin transition (^lOOMHz). Therefore, 
only transitions for which ASz — may participate. 
Considering the 9 Zeeman eigenstates, the Sz = ±2 
states are energetically isolated, the Sz = ±1 manifolds 
induce nuclear transitions but are unpolarized, and the 
Sz=0 manifold induces nuclear transitions and is polar- 
ized. The energy eigenstates with Sz=0 have energies: 

Eh. = -A + r! + r 2 (3) 
e 2 = ~(-A - V2^ 2 + (A-r 1 - r 2 )2 + v 1 + r 2 X4) 

e 3 = 1 -(-a + + (A - ri - r 2 )2 + Ti + r 2 X5) 
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Since (Fi+r2) is large compared to both the nuclear Zee- 
man energy, fiwj, and the magnitude of the dipolar inter- 
action, A, only transitions between E\ and E 3 have an 
appreciable probability of yielding the necessary energy 
for the 100 MHz nuclear transition. In the two-defect 
Zeeman basis, the state vectors for these states are: 



|1> = ^H-1,1) + 



(6) 



13} - w (| - 1,1) +«|0,0> + |1,-1» (7) 

where N 3 = V2 + a 2 and a = -2E 2 /A. Since the |0,0) 
Zeeman state is known to be selectively populated by 
optical pumping j2U the energy eigenstate |3) will become 
selectively populated. The energy eigenstates |1) and 
1 3) are taken to be a pseudo-two level system, and are 
assigned a temperature: 



T 



E 3 — Ei 
fcBln(ft) 



(8) 



For defect pairs oriented such that E 3 > E\, selective 
population of state |3) will result in a negative tempera- 
ture, and thermal contact of this pseudo-two level system 
with the nuclear spins drives the nuclear spins to a neg- 
ative temperature. After calculating transition rates, it 
is apparent these negative temperature transitions dom- 
inate, accounting for the observed NMR signal. 

The preceding model may be further quantified by con- 
sidering the rates of transitions between the pseudo-two 
level system formed by pairs of NV- centers and the 13 C 
nuclear spins. Towards that end, we first must identify 
the relationship between the known populations of the 
three Zeeman eigenstates (P\+i) , P\o) , -P|-i)) formed by 
the single NV- center, and the ratio of the two-defect 
populations, Pi/P 3 . Consideration of the 9x9 density 
matrix for coupled NV- centers yields 



Pi - P\+i)P\-i) 

p 3 = ^2 (^l+D^I-D+^l 



(9) 
(10) 



The dipole-dipole interaction between NV- centers and 
13 C spins is treated as a perturbation that induces tran- 
sitions between NV- center energy levels. Fermi's Golden 
Rule gives the direct relaxation rate for a nucleus, m, near 
a central NV- center, j, which is also coupled to another 
NV- center, k: 



W j 



km 



2tt, 



(11) 



V jm _ M7j2 7w7j ^(_| gin (0) cos (0) e -i^ jg the Hamil _ 

tonian matrix element for the transition, from Hxs- An- 
gular dependence is removed by spatial averaging 2 ^ to 
obtain \Vi m \ = —£^fi 2 jNv7i^3- In order to obtain an 
average transition rate, a sum over all neighboring NV- 
centers, k, is approximated as an integral: 



W jm = 



p{A)g(A)W 3mk {A)dA (12) 



where g is the normalized distribution of values of A, and 
p is the probability that the NV- center pair is in one of 
the states participating in the transition. For a dilute 



spin system, g(A) 



where S 



tt(<5 2 +A 2 ) ' 

and n is the number density of NV- centers^!. The inte- 
gral may then be evaluated: 



£ p(A n )Wi mk (A n )g(A n ) 



(13) 



n=l,2 



where A n are the two roots of the argument of the delta 
function in W. 

In the present model, direct relaxation of 13 C nuclei oc- 
curs only near NV- centers. Thus bulk 13 C signals build 
up due to transport of the polarization via nuclear spin 
diffusion. The governing equation for such spin transport 
is: 

^ =D ll {r 2 ^ 
dt r 2 dr dr 



where D — 6.7 * 10~ 15 cm 2 /s is the spin diffusion coeffi- 
cient due to the homonuclear dipolar interaction Hzz^, 
neglecting hyperfme "shifts" near the defectP. P is the 
nuclear spin polarization, Pi is the nuclear polarization 
in equilibrium with the temperature in Eq. (8), Wi is 
the direct relaxation rate and fi is the fraction of NV- 
center pairs in a given orientation, I. Wq is the total 
transition rate within the unpolarized Sz = ±1 mani- 
folds. Nuclear spin-lattice relaxation is neglected in the 
calculations. Using the initial condition P(t = 0) = 
for saturation of the spin transitions, and the bound- 
ary conditions ^(r = 0) = and ^(r = R) = 0, 
which follow from symmetry, a numerical solution may 
be obtained. Using n — 9 * 10 18 cm~ 3 , the numerical 
solution was rescaled and superimposed upon data taken 
at 5K in Fig. [l] This solution requires a single-defect 
NV- center polarization of more than 99% so that the 
Sz = ±1 transitions do not dominate depolarization of 
the nuclei. This polarization is higher than the 80% re- 
ported in literature^. The present sample, however, has 
a larger NV- concentration and our measurements were 
made in a larger magnetic field, such that the polariza- 
tion may be higher than reported in Ref. 21. 

The time scales apparent in Figure 1 may be compared 
to estimates assuming Eqn. 14 is governed exclusively by 
the diffusion coefficient D. The 8ppm defect concentra- 
tion presented in Ref. 18 may be used to predict the 
characteristic diffusion time, r = R 2 /D. The resulting 
time-scale, ~ 45 seconds, is of the same order of magni- 
tude as the data shown in Fig. 1. Thus, in our model, 
the contribution of NV- center dipolar fluctuations (the 
rate processes given in Eqn. 13) to the overall process re- 
flect a delicate balance of spin-diffusion, direct relaxation 
rates, and transitions within the unpolarized Sz—+/- 1 
levels. 

The equilibrium behavior as described by the spin- 
temperature treatment is consistent with the anti- 
thermal sign of the NMR signal. The polarization ki- 
netics are consistent with direct relaxation near defects 
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with spin diffusion transport of polarization to the bulk 
material. The origin of the temperature dependence of 
the ultimate polarization is not as clear. For many nu- 
clear polarization methods, temperature dependance is 
due to either the Boltzmann factor or "leakage" of po- 
larization via nuclear spin-lattice relaxation. The lack 
of temperature dependence of the polarization kinetics is 
not consistent with this picture. Instead, we propose the 
temperature dependence is due to a decreased laser pene- 
tration depth at higher temperatures. Since irradiation is 
occurring within the phonon sideband of the NV- center 
optical absorption spectrum^, an increase in absorption 
with increasing temperature is expected as more phonons 
are available to participate in the absorption, leading to a 
decrease in optical penetration and a smaller irradiation 
volume. Thus, a model for the temperature and power 
dependence, while beyond the scope of this article, could 
consider their effect on the spatial profile of light inten- 
sity I(z), and in turn, the polarization. 

The present model, as summarized in Equation 14, 
does not account for all contributions to the spectral 
density of fluctuating fields at the 13 C nucleus. The 
temperature-dependence of the 13 C lineshapes suggests 
that the are other factors influencing the spin dynamics 
in this sample, likely other paramagnetic defects such as 
substitutional nitrogen, which have highly temperature 



dependent polarization and are also coupled to the nu- 
clei and may contribute to relaxation processes. This and 
other effects, for example, the dynamics of NV- center 
photophysics, would effectively result in a modification of 
the function g(A). However, our simple model accounts 
for the dominant features of the data, and may be further 
tested by noting that the temperature in Eq. Q depends 
on the crystal orientation through the parameter T. 

We consider the results of this study to be important to 
several scientific communities. The mechanism described 
herein portends control over the nuclear spin bath in di- 
amond, whose interaction with NV- centers is an area of 
active research. This work also suggests that bulk 13 C 
NMR measurements probe the dynamics of an ensemble 
of NV- centers, providing a methodology complimentary 
to optical (single NV- center) techniques. Spin-polarized 
nuclei in diamond may also be useful for sensitivity- 
enhanced NMR techniques where polarization localized 
near the crystal surface may be transfered to other nuclei 
(e.g. adsorbed molecules) for enhanced NMR sensitivity. 
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